MKPH-T-93-04, IU/NTC 92-20 March 24, 1993 



Nuclear Shadowing in a Parton Recombination Model 

S. KUMANO * 

Institut fur Kernphysik, Universitdt Mainz 
6500 Mainz, Germany 
and 

Nuclear Theory Center, Indiana University 
Bloomington, Indiana 4^408, U.S.A. 



ABSTRACT 

Deep inelastic structure functions F% (x) are investigated in a Q 2 rescaling model with 
parton recombination effects. We find that the model can explain experimentally measured 
F^x) structure functions reasonably well in the wide Bjorken— x range (0.005 < x < 0.8). 
In the very small x region (x < 0.02), recombination results are very sensitive to input 
sea-quark and gluon distributions. 
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1. Introduction 

Several years have passed since the discovery of nuclear modifications in structure 
functions ^(x) by the European Muon Collaboration ("old" EMC effect) [1]. In spite 
of initial expectation of an explicit quark signature in nuclear physics, it is still not clear 
whether the effect originates in nucleon substructures or just in nuclear physics. There are 
several proposed models for explaining the EMC effect. Some models tried to interpret it 
by explicit quark effects, e.g. Q 2 rescaling models [2] and six-quark bag models. Others 
investigated models based on conventional nuclear physics, such as nuclear binding and 
pion excess. For details of these models, we suggest that interested readers look at summary 
papers in Ref. 3. Although there are still experimental activities [4] for investigating these 
different models, it is rather difficult to discriminate among these models. 

Most of these investigations discuss a "global" EMC effect in the sense that the effect 
is averaged over all constituents in a nucleus. However, it is interesting to investigate 
possible semi-inclusive or semi-exclusive processes for finding a "local" EMC effect [5] and 
possible relations between a local gluonic EMC effect and the J/ip suppression [6]. 

Considering the fact that the average nucleon separation in nuclei is 2.2 fm and the 
nucleon diameter is 1.8 fm, we expect that nucleons in nuclei could overlap strongly. If 
a multiquark cluster is formed in a nucleus due to this overlap, the confinement radius 
for such a quark should be different from the one in a free nucleon. Using this kind 
of simple picture in a small Q 2 region and the Q 2 -evolution equation to compare with 
experimental data at large Q 2 , we obtain a simple prescription for the nuclear structure 
function F^(x,Q 2 ) [2]. It is related to the nucleon structure function by a simple Q 2 
rescaling, F^(x, Q 2 ) — F^(x, i A Q 2 ), where £ A is called as rescaling parameter. Although 
the simple size change could be too simple to explain many details of nuclear physics [7], 
it is a useful effective model in explaining deep inelastic data in the medium x region. 

The New Muon Collaboration (NMC) [8] and the Fermilab E665 collaboration [9] 
recently measured accurately the structure functions ^(x) at very small x. These data 
as well as EMC data [10] provide an opportunity for investigating physics details in the 
shadowing region (x < 0.1). The shadowing means that central constituents are shielded 
due to surface constituents, hence the cross section behaves like A a (a < 1). There are 
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different ideas for explaining the shadowing phenomena: vector meson dominance [11], 
parton recombination [12—17], qq fluctuations of the virtual 7 [18], pomeron dominance 
[19] (with Pauli blocking [20]), and hybrid models [21]. 

At first, there are models based on the traditional idea, the vector meson dominance 
model [11]. A virtual photon transforms into vector meson states (p, u, 0), which then 
interacts with a target nucleus. The propagation length of the hadronic (H) fluctuations 
is given by A « 1/\Ejj — E^\ = 2i//(M^ + Q 2 ) 0.2/x fm. For x < 0.1, the propagation 
length (>2 fm) exceeds the average nucleon separation (2.2 fm) in nuclei and the shadowing 
takes place due to multiple scatterings. For example, the vector meson interacts elastically 
with a surface nucleon and then interacts inelastically with a central nucleon. Because this 
amplitude is opposite in phase to a one-step amplitude for an inelastic interaction with 
the central nucleon, the nucleon sees a reduced hadronic flux (namely the shadowing). 

The parton recombination model [12—17] has been investigated as a mechanism for 
explaining the shadowing within the framework of a quark-parton model. In an infinite 
momentum frame, the average longitudinal nucleon separation in a Lorentz contracted 
nucleus is L = (2.2 j "m)M 'a/ 'F 'a = (2-2 fm)m N /p N: and the longitudinal localization 
size of a parton with momentum xp N is AL = l/(xp N ). If the parton dimension (AL) 
exceeds the average nucleon separation (L), partons from different nucleons could interact 
(we call this parton recombination or parton fusion) significantly and the shadowing could 
occur due to processes in Figs, lb, lc, le, and If. Partons with momentum fraction x are 
lost in these processes, so that their contributions are negative AF2(x) < (shadowing). 
However, Ai^x) depends much on input sea-quark and gluon distributions as we find in 
sections 3 and 4. A relevant x region for the shadowing is obtained by using L < AL and 
we find x < 0.1. Even though the recombinations could produce shadowing type effects in 
the small x region, they are very small compared with experimental shadowings if they are 
calculated at Q 2 =b GeV 2 . On the other hand, modifications of ^(x) due to recombination 
effects on gluon distributions are large [17]. Details of numerical results are discussed in 
section 3. 

Because the recombination results by Close, Qiu, and Roberts [17] are not compared 
in detail with the recent experimental data, we investigate whether the model can explain 
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the EMC, NMC, and E665 shadowing data in this investigation. Furthermore, we study 
whether the Q 2 rescaling model combined with parton recombination effects [14,15] can 
explain the F 2 (x) structure functions in the whole Bjorken— x range [16]. 

Earlier comparisons with experimental data have been made in similar parton pictures 
[15]. Let us clarify differences from and advantages over the previous investigations. (1) 
Covolan, Predazzi, and others [15] used the Qiu's parametrization [13] for the shadowing 
and combined its results with the rescaling. We note that the Q 2 evolution is well 
investigated in Ref. 13; however, the x dependence is rather assumed. Therefore, we cannot 
have a dynamically consistent picture by using the shadowing-parametrization with the 
rescaling. The x dependence of the shadowing was later investigated in Ref. 17. Using this 
shadowing picture combined with the rescaling, we could possibly obtain a dynamically- 
consistent parton model in the wide x-range. This is the purpose of our investigation. 
(2) We study Q 2 evolution effects in the combined model. (3) We investigate effects of 
input sea-quark and gluon distributions on the shadowing. (4) Improved experimental 
data became available by the NMC and E665. 



If a nucleon is in a nucleus, parton distributions are modified because of the existence 
of neighboring nucleons. Partons from different nucleons could interact with each other, 
and the interactions become important, especially in the shadowing region. Processes 
contribute to modifications in a quark distribution q(x) are shown in Fig. 1. For example, 
Fig. la indicates that a quark from the nucleon 1 fuses with a gluon from the nucleon 2 
and produces a quark with momentum x. Because the process creates a quark with the 
momentum x, this is a positive contribution to the quark distribution at x. There are five 
other contributions as shown in Fig. 1. In general, a modification of a parton distribution 
Psixs), due to the process of producing the parton p% with the momentum X3 by a fusion 
of partons p\ and p 2 , is given by [17] 



Af>3(>3) = K J dx\dx 2 p\{x\) p 2 (x 2 ) T pip2 ^ P3 (xi, x 2 , x 3 = xi + x 2 ) S(x, xi, x 2 ) , (1) 
where K is given by K = 9A 1 / 3 a s /(2R 2 l Q 2 ) [22]. Nuclear radius R is R = RqA 1 / 3 with 
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i?o=l-l frn [23] and the strong-interaction coupling constant is a s (Q 2 ) = 4n /[/3oln(Q 2 / A 2 )] 
with fio = 11 — 2iVj/3, Nf=3. The momentum conserving 5 function is given by 
5(x — x\ — X2) for Figs, la and Id and 5{x — x\) for Figs, lb, lc, le, and If. The parton 
fusion function T{x\, X2, X3) is a probability for producing a parton p% with momentum 
£3 by a fusion of partons p\ and P2 with momenta x\ and xi respectively. Four possible 
parton fusion processes are shown in Fig. 2. It is related to a splitting function P piP3 (z) 
in the Altarelli-Parisi equation [24] by 

£3 £3 

The splitting function P PlPs (z) is the probability of finding the parton p\ with fraction z 
of the parent momentum in the parton ^3. C PlP2 ^ P3 is the ratio of color factors in pro- 
cesses p!p 2 -> P3 and ps -> pip2- For example, C gG ^ g = Yl (tfkYtfk = 

(l,a),k (k),l,a 

[Cf/(Nq — 1)]/Cfi where averages are taken over initial color indices (k in the denomina- 
tor, / and a in the numerator). tf k is given by the SU(3) C Gell-Mann matrix by tf k = Xf k /2, 
and C F is C F = (N 2 - 1)/(2N C ) [25]. In this way, the color factors are calculated as: 

r - 1 1 r _ _ 8 

qG ^ q ~ N 2 - 1 ~ 8 ' qq—>G - - Q , 

where Tp = 1/2 and iVc = 3. Using the above results and splitting functions in Ref. 24, 
we obtain the parton fusion functions: 

^qG^q'(x 1 ,X 2 ,X 3 ) = ^Z(l+Z 2 ) , 

r q q-+G(xi, X 2 , X 3 ) = ^ Z (1 - 2) [Z 2 + (1 - ^) 2 ] , 

r G ^(xi,x 2 ,x 3 ) = ^ (1-2;) [l + (l-z) 2 ] , (4) 

^G 1 G 2 ^G 3 {X1,X2,X 3 ) = \ Z {I- Z) [—?—+- - + Z(1-Z)} , 

4 1 — Z Z 

where x$ = x\ + X2 and z = x\/x%. 
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Using these fusion functions, we calculate modifications of quark distributions in nuclei 
for processes shown in Fig. 1. Explicit expressions for the modifications Aqi(x) are shown 
in Appendix (Eq. (Al)). From Eq. (Al) and a similar expression for Aq~i(x), we obtain 
recombination contributions to the structure function F2(x) in Eqs. (A2.1) — (A2.4). Direct 
recombination effects on the EMC effect are given by the ratio 

AF A (x Q 2 ) 

^(recombination) = 1 + 2 ' , (5) 

r 2 [ x i Q ) 

where AF<f(x, Q 2 ) is calculated by using Eqs. (A2.1) — (A2.4) with input parton distribu- 
tions at Q 2 . The most important physics for the shadowing in the recombination model is 
the effect due to modifications in gluon distributions. Recombination effects on the gluon 
distribution are given in Eqs. (A3. 1) — (A3. 4). Calculating these equations, we find that 
the recombination mechanism produces large shadowing effects on the gluon distribution 
[17]. 



3. Gluon-shadowing effects on F^(x) 

3.1 An approximate way to take into account the gluon shadowing 

Gluons and quarks are coupled to each other and their distributions are related by 
the Altarelli-Parisi equation. Because of the large modifications in G(x < 0.1) due to 
the recombinations, we expect that sea-quark distributions are also affected by the gluon 
modifications. In the very small x region, the Altarelli-Parisi equation is dominated by 

gluon dynamics, and it can be solved analytically. As a result, there is a relationship 

x d 

between a sea quark and gluon distributions [26]: xqf ea (x) = — — — [xG(x)], where z=u, 
d, or s. Using this relation, we obtain a contribution to F2(x) from the modification in the 
gluon distribution: 

5F 2 (x, Q 2 ) = -0(x o - x) 1 ^i xAG ( x > Q 2 )] . ( 6 ) 

where xAG(x, Q 2 ) is given in Eqs. (A3. 1) — (A3. 4), and a step function 9(xq — x) is 
introduced because the relation is valid only at very small x. It is defined by 9(xq — x) = 1 
(for x < xq) or (for x > xq). We should note that this is a crude estimate of the gluonic 
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contribution. To be precise, the recombination calculations are done at small Q , then 
the results are evolved to the Q 2 region where experiments were done, by using a nuclear- 
modified Altarelli-Parisi equation [13]. This approach is discussed in section 4. Combining 
all the contributions, we obtain nuclear structure functions F 2 A (x) calculated in our parton 
model (the Q 2 rescaling model with parton recombination effects): 



where F A is given by F A (x, Q 2 ) — F 2 (x, i A Q 2 ) in the rescaling model. 

3.2 Results by using the analytical solution at small x 

We calculate nuclear structure functions and investigate whether our theoretical results 
are compatible with experimental data [27] by a SLAC group [28], EMC [10,29], NMC [8], 
and E665 [9] . It is interesting to investigate not only the shadowing region at small x but 
also the large x region, which is traditionally called "nucleon- Fermi- motion" part. 

In evaluating equations in Appendix, we first assume that a leak-out parton is a 
sea quark or a gluon and that the momentum cutoff function [17,30] for this parton is 
taken as w(x) = exp(—m 2 N ZQX 2 /2), namely q*(x) = w(x)q sea (x), q *(x) = w(x)q(x), and 
G*(x) = w(x)G(x). Input parton distributions are given by a recent parametrization by 
Martin, Roberts, and Stirling (MRS) [31] or by Kwiecinski, Martin, Roberts, and Stirling 
(KMRS) [32]. Q 2 = 5 (or 4) GeV 2 is used in the input parton distributions and for 
calculating K in Eq. (1). The cutoff function w(x) supposedly takes into account effects 
of bound partons and it is shown as the function of x in Fig. 3. As discussed in the 
introduction, x = 0.1 corresponds to the length scale 2 fm, the average nucleon separation 
in nuclei. Hence, we expect that only partons p* with x^O.l could participate in the 
recombinations. We then find in Fig. 3 that an appropriate choice is zq=2— 3 fm. We 
show theoretical results with zq=2 fm in this section; some zq dependencies are discussed 
in section 4.3. xq in Eq. (6) is taken as x$ = 0.1. We do not show the ratio in Eq. (7) in 
the x range 0.1 < x < 0.2 because physics is not well described by the crude estimate in 
Eq. (6) for taking into account the gluonic modifications. This x region is investigated by 



Ff{x, Q 2 ) 
F 2 D (x,Q 2 ) 



F A (x, Q 2 ) + AF 2 A (x, Q 2 ) + 5F A (x, Q 2 ) 
F 2 D (x,Q 2 ) 



(7) 
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a more-complete Q 2 -evolution picture in section 4. The QCD scale parameter A in a s (Q 2 ) 
is taken as A=0.2 GeV. 

In our theoretical analysis, Ag, Sn, and Xe targets are assumed as 107 Ag, 118 Sn, and 
131 Xe nuclei. We take the rescaling parameters in Table II of the Close, Jaffe, Roberts, and 
Ross paper [2]. They are £ A =1.43 ( 4 He), 1.60 ( 12 C), 1.86 ( 40 Ca), 2.17 ( 107 Ag), and 2.24 
( 118 Sn). £ A =2.24 is assumed for the 131 Xe nucleus. The above parameters are obtained at 
Q 2 =20 GeV 2 . Even though the model can explain the data, it is known that the original 
fits are no longer valid if we include Fermi-motion effects [33] . Furthermore, the rescaling 
parameters are calculated in a "semi-classical" way in the sense that overlapping volumes 
are estimated simply by the geometrical ones. Considering these points, we do not think 
that the parameters are not well defined in their overall magnitudes. According to the 
dynamical rescaling model by Close, Roberts, and Ross [2] in 1988, nuclear dependence is 
the result of a scale change. So, we could in principle set a rescaling parameter by fitting 
F 2 (x) data of a medium-size nucleus, and then use the density dependence of the rescaling 
model in calculating £ A for other nuclei. However, we find in our numerical analysis that 
the above £ A explain the SLAC data reasonably well, so that we simply keep the original 
rescaling parameters in Table II of Ref. 2. 

We discuss results for the calcium nucleus in Figs. 4a and 4b. In Fig. 4a, theoretical 
results in Eqs. (5) and (7) are compared with the SLAC data [28] in the linear x scale. For 
simplicity, the deuteron structure functions are assumed as F 2 D (x) = {F%(x) + F£(x)]/2 = 
F 2 (x) in our theoretical analysis. Two sets of results are shown: the solid curves are 
obtained by using the MRS-1 distributions [31] and the dotted curves are by the KMRS- 
B0 [32]. A, B, C, and C in the figure show 

(A) [F 2 N (x, Q 2 ) + AF 2 (x, Q 2 )]/F 2 N (x, Q 2 ) , 

(B) [F 2 N (x, £ A Q 2 ) + AF 2 (x, Q 2 )]/F 2 N (x, Q 2 ) , 

(C) [F 2 N (x, Q 2 ) + AF 2 (x, Q 2 ) + SF 2 (x, Q 2 )]/F 2 N (x, Q 2 ) , 

(C) [F 2 N (x, ^ A Q 2 ) + AF 2 (x, Q 2 ) + 5F 2 (x, Q 2 )]/F 2 N (x, Q 2 ) , 

where AF 2 (x,Q 2 ) and 5F 2 (x : Q 2 ) are recombination results with Q 2 -rescaled input dis- 
tributions. Direct recombination contributions, shown by the curves with A, are positive 
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in the medium and large x regions. These come from the quark-gluon fusion processes in 
Figs, la and Id. Combined results with the Q 2 rescaling are shown by the curves with B. 
The C and C curves are explained in the following paragraph. Fig. 4a indicates that our 
model results are in overall agreement with the experimental data. The EMC effect in the 
medium x region is mainly due to the Q 2 rescaling in our model, because the recombination 
contributions are rather small. It is noteworthy in these figures that the recombination 
can explain nuclear structure functions in the large x region without explicit Fermi-motion 
effects. Physics in this region is usually attributed to the nucleon Fermi motion in the 
nucleus. In our model, a quark with x^l could be pushed by a gluon from a different 
nucleon and becomes a quark with x > 1. This is the reason why the ratio goes to infinity 
at x — > 1 in Fig. 4a. Perhaps, there exists a correlation between these two apparently 
different descriptions: the parton recombination and the nucleon Fermi motion. 

In Fig. 4b, theoretical results are compared with EMC [10] and NMC [8] data in 
the logarithmic x scale. In the small x region, the processes in Figs, lb, lc, le, and If 
significantly contribute and produce negative (shadowing) results in the case of KMRS-BO. 
However, they are too small to explain the whole shadowing data if the recombinations are 
calculated at Q 2 =h GeV 2 . The important point in our model is the shadowing which is 
produced through the gluon distributions. Combined results are shown by the curves with 
C or C' in Fig. 4b. The figure indicates that the shadowing produced through modifications 
in gluon distributions can explain the EMC and NMC shadowing data fairly well, even 
though the direct recombination effects on i^x) (shown by A) are very small. The curves 
C are obtained without rescaling and C' is with the Q 2 -rescaled inputs, where the same 
£ A is used in all parton distributions. We find that the differences between the C and C' 
curves are rather small considering the fact that the Q 2 rescaling produces large sea-quark 
modifications. This is because the Q 2 -rescaled sea-quark distributions and contributions 
(Eq. (6)) from <2 2 -rescaled gluon distributions almost cancel each other. 

We have shown that our model can explain the NMC shadowing in Fig. 4b. However, 
the results are very sensitive to input gluon distributions in the small x region. Because 
the E665 collaboration measured F* e {x) / 'F^ \x) at very small x recently, we compare our 
results with the E665 data. In order to find details of input gluon distribution effects, we 
use parton distributions with "hard", "soft", and u l/y/x n gluon distributions given by the 
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MRS [31]. Analytical distributions at Q 2 =4 GeV 2 are used as the input distributions. The 
MRS soft, hard, and 1/y/x gluon distributions are given by 5 G = 0, r\ G = 5, 7 G = 0, 
5 G = 0, r\ G = 4, 7 G = 9, and 5 G = —1/2, r\ G = 4, r y G = 9 respectively for the 
notation xG(x) = Aqx s g (1 — x) v g(1 + 7 G x). These gluon distributions are shown in 
Fig. 5a and they are very different especially at x < 0.01. Our shadowing results are 
shown in Fig. 5b. Because the experimental Q 2 of the E665 data becomes very small 
(Q 2 « 1 GeV 2 ) at x < 0.001, our perturbative results cannot be compared with the 
E665 data. This is the reason why curves in the region x < 0.001 are shown by the dashed 
ones. As shown in the figure, our calculations produce widely-different shadowing results 
depending on the input gluon distribution. If it is the 1 / \fx distribution, shadowings are 
too large and are contradictory to the experimental data. If it is the hard one, theoretical 
results underestimate the shadowing. Our theoretical results are nearly consistent with 
experimental data if the input distributions are the soft MRS distributions at Q 2 —A GeV 2 . 

We found that our model with the recombination and the rescaling works in the wide- 
x region, if an appropriate input gluon distribution is taken. Nevertheless, it is not a 
complete description in the sense that F^(x) in the x ~ 0.1 region cannot be calculated. 
Furthermore, a curve in the x > 0.2 region and one in the x < 0.1 are not smoothly 
connected in Fig. 4b. It implies that the analytical solution Eq. (6) at a; < 0.1 is used 
beyond the applicable region. In order to have detailed comparisons with experimental 
data and to have a dynamically consistent picture, Q 2 evolution effects should be fully 
investigated. We address this point in the following section. 
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4. Q 2 evolution of nuclear structure functions 

4.1 Q 2 evolution 

The approximate description with Eq. (6) is good enough for rough magnitude 
estimates. For detailed comparisons with the experimental data in the x w 0.1 region, 
we should solve (Revolution equations. Namely, we calculate the recombination and the 
rescaling at small Q 2 (= Qq) : then obtained distributions are evolved to those at larger 
Q 2 where experimental data are taken. However, this is not an easy work by the following 
reasons. 

First, the shadowing data are mostly taken in the Q 2 region of a few GeV 2 . We need 
to assign a very small value to Q 2 , for example Q 2 =l GeV 2 , in order to study evolution 
effects by comparing our results with the existing data. The perturbative QCD, especially 
in the leading order, would become questionable in such small Q 2 region. Nevertheless, we 
present the evolution effects in this report because the gluon shadowing connected with the 
Q 2 evolution is the essential part of our model in explaining the F<i shadowing, as we found 
in the previous section. Next-to-leading (and higher) order corrections are important in 
the small Q 2 and small x region. Hence, our investigation should be considered as a first 
step for understanding F^(x) dynamically in the wide x region within the framework of a 
quark-parton model. 

Secondly, the evolution equation for parton distributions in a nucleus is given in 
Ref. 13. Finding a numerical solution for the equation is by itself a significant research 
problem. So, instead of solving the nuclear evolution equation, we use the Altarelli-Parisi 
equation for our Q 2 evolution. Both evolution results are qualitatively similar; however, the 
nuclear-evolution shows smaller <3 2 -dependence for large nuclei [13]. Because our results 
indicate small Q 2 dependence for large nuclei in section 4.2, this problem is not considered 
to be very serious. 

We calculate F^(x) in the following way. Valence-quark distributions in the nucleon 
are modified according to the Q 2 rescaling mechanism at Q 2 : V{x, Qq) = V(x, £,\Qq)- 
In order to satisfy the momentum conservation J dxx[u v (x) + d v (x) + S(x) + G(x)] = 1, 

sea-quark and gluon distributions are simply modified as S(x, Q 2 ) = C sg S(x, Qq) and 
G{x, Q 2 ) = C sg G{x, Qq). (It is not clear whether the rescaling picture could be used 
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also for sea-quarks and gluons. This issue is discussed in section 4.3.) Obtained parton 
distributions are then used for calculating the recombination effects in Appendix, and we 
get nuclear parton distributions at Qq. These distributions are evolved to those at larger 
Q 2 by using the ordinary Altarelli-Parisi equation [24]. For the flavor-singlet case, it is 

l 

j-qs&t) = / [q s (y,t) P qq (-) + G(y,t) P qG {-) ] , 
j y y y 

X 

1 

§-G(x,t) = [ ^ [ q s (y,t) P G A + G(y,t) P GG £) ] , (8) 
vi* j y y y 

X 

where q s is the flavor-singlet distribution given by q s = ^^(qi + <?i), Pij is the splitting 

i 

function, and t is defined by t = — (2/ f3o)ln[a s (Q 2 ) /a s (Q 2 )]. We used a solution of this 
integrodifferential equation in the leading order [34]. In the flavor-nonsinglet case, the 
gluon does not enter into the evolution equation: 

l 

— q ns (x,t) = — q ns (y,t)P qq (-) , (9) 
ot J y y 

X 

where q ns = qi — qi. We use the distributions at Qq calculated by the rescaling 
and the recombinations in the above equations: q ns (x,Q^) = u v (x,Qq) + d v (x,Q^), 
q s (x,Q 2 ) = u v (x,Q 2 ) + d v (x,Q 2 ) + S(x,Q 2 ), and G(x, Q 2 ) = G(x,Q 2 ). We then 
obtain evolved distributions by solving the above equations, and F2 at Q 2 is given by 
F2{x,Q 2 ) = [xq ns (x, Q 2 ) + 4xq s (x, Q 2 )]/18. The use of the above equations instead 
of the nuclear evolution equations causes a problem in the large x region, because the 
recombinations produce distributions at a; > 1. However, it is not a serious problem in 
comparison with the experimental data discussed in the next section. 

4.2 Q 2 evolution results 

Structure functions F<f(x) are calculated in the following way. Because the parton 
distributions in the MRS have simpler analytical forms than those in the KMRS, we first 
employ the MRS-1 (soft gluon) as our parton distributions in the nucleon. Using the 
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analytical form of the MRS-1 at Q 2 =A GeV 2 and the evolution subroutine in Ref. 34, 
we calculate distributions at Q 2 (~ 1 GeV 2 ). Then, the Q 2 rescaling is used for valence 
distributions and the constant C sg is determined so that the momentum conservation 
is satisfied. Input distributions at Q 2 are parametrized by a simple analytical form, 
xp(x) = Ax s (l — + 7a;). Constants A, S, r], and 7 are obtained by fitting numerical 

results. The recombinations are calculated by using the modified input distributions at 
Q 2 , and obtained results are again fitted by xp(x) = Ax 5 (l — x) v (l +7^). For the valence 
quarks, we used xV(x) = Ax 5 (l — /3x) v (l + 72;), where j3 < 1, so that we could take into 
account the fusion effects at large x (=0.8—0.9). The recombinations produce distributions 
at x > 1; however, we simply neglected such distributions in our Q 2 evolution. This causes 
slight inconsistency, but the neglected effect on the momentum conservation is of the order 
of 10~ 4 %, which is insignificant. 

We discuss results for the nuclei He, C, Ca, Ag, Sn, and Xe in Figs. 6-9. In Figs. 
6a-9a, our theoretical results are compared with SLAC [28] and EMC [29] (E665 [9]) data 
in the linear x scale. In order to see the shadowing due to the Q 2 evolution, we should 
choose Q 2 ~1 GeV 2 because the data were taken in the region of a few GeV 2 . The value of 
Qq 1 is adjusted so that our results agree with the shadowing data for 40 Ca, and we obtain 

2 r, 

<5o=0.8 GeV . We fixed the cutoff at zo=2 fm as it was used in section 3. We discuss 
dependence of our results on these parameters in the next section. Another problem is 
how we should choose the rescaling parameters at very small Q 2 . As we did in section 3, 
we may simply determine £ A for a medium-size nucleus and then use the A-dependence for 
other nuclei. However, it turns out that an original rescaling parameter for a medium-size 
nucleus is not a bad choice, so we simply use the parameters in Ref. 2. Three curves with 
<3 2 =0.8, 5.0, 20.0 GeV 2 are shown in Figs. 6-9. Figures 6a-9a indicate that our model 
can explain the experimental data well except for the x ~ 0.7 region of large nuclei. The 
results in this region depend much on the input valence distribution and the cutoff (see 
section 4.3). If we change the input to the KMRS-B0 distributions, our calculations work 
very well for large nuclei, but the EMC effect for 4 He is overestimated. 

In Figs. 6b-9b, theoretical results are compared with EMC [10], NMC [8], and E665 
[9] data in the logarithmic x scale. Our results are not shown at x < 0.001 in Fig. 9b 
because the used evolution subroutine [34] does not have good accuracy in such region and 
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our perturbative calculations cannot be compared with the data with very small Q 2 . We 
find in Figs. 6b-9b that our results agree quite well with the experimental data. In the 
small nuclei (Figs. 6b and 7b), there are some Q 2 dependence at Q 2 = 0.8 — 2 GeV 2 , but 
it becomes Q 2 independent at Q 2 > a few GeV 2 . We have very small Q 2 dependence in 
the larger nuclei as shown in Figs. 8b and 9b. We note that rather large Q 2 dependence 
in Figs. 6b and 7b should not be taken seriously at this stage, because small variations in 
input sea-quark and gluon distributions change the dependence drastically as we find in 
section 4.3. 

4.3 Dependence on the parameters 

We discuss how our results depend on the parameters, the input distributions, and the 
rescaling assumption. We first check the choice of Q 2 . We choose Q^=2.0 GeV 2 instead 
of 0.8 GeV 2 in Fig. 10. Distributions at 2 GeV 2 are calculated and are then evolved to 
those at 5 and 20 GeV 2 . It is obvious that our calculations underestimate the experimental 
shadowings and that the Q 2 dependence is very different from the one in Fig. 8b. The 
differences are due to the Q 2 factor in K (see Eq. (1)) and due to input sea-quark and 
gluon distributions. In fact, we find in the following that slight modifications in those 
distributions change our shadowing results at Q 2 significantly. 

In order to check input-distribution effects, we take the KMRS-BO distributions which 
are slightly different from those of the MRS-1. These distributions are shown in Fig. 11. 
Even though two distributions do not look very different, Q 2 dependency of the KMRS-BO 
results for the carbon nucleus in Fig. 12a is quite different from that in Fig. 7b. The 
KMRS-BO results for the calcium in Fig. 12b are rather similar to those of the MRS-1 in 
Fig. 8b. Because the parton distributions in the small- x region are not well known, we 
inevitably have such uncertainty in our model. 

Effects of a momentum cutoff for leak-out partons are shown in Figs. 13a and 13b. 
From Fig. 3, we decided to take zq=2 fm and we have been using this cutoff so far. Results 
with -2o=3 fm are shown by the dashed curves. Our model with larger zo produces smaller 
recombination effects at large x and larger effects at small x. The agreement with the 
data becomes better in Fig. 13a in the medium- £ region; however, it becomes worse in 
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the small-x region in Fig. 13b. Because the cutoff is one of the important factors in our 
model, efforts should be made to estimate it in an independent way. 

The last issue is the rescaling. The original rescaling model is intended for valence 
quarks. It is not clear whether sea-quarks and gluon follow the same rule of the scale 
change in nuclei. In our analysis of this section, we used the rescaling only for valence 
quarks. Sea-quark and gluon distributions are simply increased by a constant amount 
iC sg ) so that they carry a momentum deficit produced by the valence rescaling. This is 
an arbitrary assumption. We may choose other scenario, for example the Q 2 rescaling for 
all partons. Using this rescaling picture, we obtain the results in Fig. 14. The shadowings 
are much underestimated in this case. Although numerical values do not agree completely 
with those in Fig. 4b, we find a similar tendency in the sense that calculated shadowings 
become smaller in the x region (0.01 < x < 0.1) if the rescaling is used for all partons. 

We learned the following from the above analyses. There are a few parameters and 
assumptions in our model and obtained results depend inevitably on their choice. The 
major factors are (1) Q 2 , (2) input S(x) and G(x) in the small x, (3) zq : and (4) Q 2 
rescaling for (valence) partons. We find that magnitude of our shadowing and the Q 2 
dependence are very sensitive to the above factors. Nevertheless, we find that it is possible 
to choose a set of the factors within our initial-expectation ranges (e.g. Q 2 ~1 GeV 2 , 
zq ~ 2 fm) and we explain many existing data. Obviously, much efforts should be done 
for investigating these factors and also next-to-leading order effects. 

4. Conclusions 

We investigated nuclear structure functions i 7 ^, Q 2 ) from small x(m 10~ 3 ) to large 
x(~ 0.9) in a parton model based on a Q 2 rescaling model with parton recombination 
effects in order to compare them with the recent experimental data. As a result, we 
obtained reasonably good agreement with the experimental data in the region (0.005 < 
x < 0.8). In the large x region, the ratio (F^ (x) / (x) > 1) is explained by quark-gluon 
recombinations, which produce results similar to those by the nucleon Fermi motion. In 
the medium x region, the EMC effect is mainly due to the Q 2 rescaling mechanism in 
our model. In the small x region, shadowing effects are obtained through modifications in 
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gluon distributions. However, our shadowings at very small x(< 0.02) are very sensitive to 
the input gluon distribution. We have a few parameters in our model; however, we could 
choose a set of the parameters and explain many existing experimental data. Q 2 variations 
of our shadowing results depend much on the input sea-quark and gluon distributions, and 
also on the parameters. 
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Appendix 

Although detailed formalisms are given in Ref. 17, explicit equations used in this 
investigation are shown in the following. Using Eqs. (1) — (4) and changing the integration 
variable x\ <-> X2 for the process in Fig. le, we obtain modifications of quark distributions 
in nuclei due to the parton recombination mechanism: 



K f dx2 
x ■ Aqi(x) = +— / X2G*(x2) 



6 j 12 





(x - x 2 )qi(x - x 2 ) {l + ( 
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(x + x 2 y 



(Al) 



where * indicates a leak-out parton. For example a gluon, which leaks out from the nucleon 
2, interacts with a quark in the nucleon 1 and produces a quark with momentum x in the 
first term of Eq. (Al). Explicit Q 2 dependence in parton distributions is not shown in 
this appendix in order to simplify the notations. Modifications in antiquark distributions 
are obtained in the similar way by changing qi <-> q\ in Eq. (Al). Using the above 
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expression, we obtain parton recombination effects on the structure functions F 2 (x) by 
AF 2 (x) = J2 e l x i A( li( x ) + Aqi(x)} 



AF 2 (x) = AF 2 [1 \x) + AF 2 {2 \x) + AF 2 {3 \x) 
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where indicates that the summation is averaged over partons in the proton and the 
neutron. For example, they are efqi = [-^(u + d) + ^«] proton and e 2 ^ = 
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5 2 
I" — (uu + dd) H ssl , . 

L]_g v 18 ■'proton 

In the similar way, modifications in a gluon distribution in a nucleus are obtained as 

[35] 

xAG(x) = xAG (1 \x) + xAG (2 \x) + xAG {3) (x) (ASA) 
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where ^ = ^(<?i) pr0 ton and ^ qiq\ = ^{^0 proton- Using these expressions, we find 
«- J _ conservation 'is e„ sa«e d , / . {A . W + 
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+xAG(x)~\ = 0. For numerical analysis, the last integral in AG^(x) should be evaluated 

1 X l 

by separating the integral region J dx<i = J dx<i + J dx2 and by changing variables 

a; 

x 

x' = x — X2 in the integral J dx2- We use the following equations for evaluating AG^\x) 



o 

in order to cancel out infinities: 
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Nuclear gluon distributions in our model are discussed in Ref. 36 and they are 
compared with recent NMC measurements of Gs n {x)/Gc{x). 
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Figure Captions 



1. Schematic pictures of parton recombination processes. 

2. Parton fusions for (a) qG — > q, (b) qq — > G, (c) Gq — > (/, and (d) GG — > G. 

3. Momentum cutoff for leak-out partons, w(x) = exp(—m 2 N z 2 l x 2 /2). 

4. Comparisons with (a) SLAC data [28], (b) EMC-90 [10] and NMC [8] data for 
Ca. Solid (dotted) curves are obtained by using the MRS-1 (KMRS-BO) input 
distributions. Q 2 =5 GeV 2 and zo=2 fm. (A) recombinations, (B) recombinations 
+ rescaling (£ A = 1.86). (C) recombinations with gluon-shadowing effects, (C) the 
same as C except for the rescaling. The only KMRS-BO curve is shown in C'. See 
text for detailed explanations of A, B, C, and C. 

5. (a) MRS gluon distributions. See Ref. 31 and text for details of hard, soft, and 
1 / \fx gluon distributions, (b) Comparisons with E665 data for Xe [9] and EMC-88 
data for Sn [29]. Theoretical results are obtained (b) for the Xe nucleus (£ A = 2.24). 
Q 2 =4 GeV 2 and z =2 fm. 

6. Comparisons with (a) SLAC data [28], (b) NMC data [8] for He. Q 2 =0.8, 5, and 20 
GeV 2 for the dotted, solid, and dashed curves respectively. zo=2 fm, £^ = 1.43 and 
Q 2 =0.8 GeV 2 . 

7. Comparisons with (a) SLAC [28] and EMC-88 [29] data, (b) EMC-90 [10] and NMC 
[8] data for C. Notations for the dotted, dashed, and solid curves are the same in 
Fig. 6. Cl = 1-60. 

8. Comparisons with (a) SLAC data [28], (b) EMC-90 [10] and NMC [8] data for Ca. 
£ = 1-86. 

9. Comparisons with (a) SLAC data for Ag [28], EMC-88 data for Sn [29], and some 
of E665 data for Xe [9], (b) E665 data for Xe [9] and EMC-88 data for Sn [29]. 
Calculated results are for Ag = 2.17) in (a) and for Xe = 2.24) in (b). 

10. Dependence on Q 2 . Q 2 =2.0 GeV 2 is taken. Q 2 =2, 5, and 20 GeV 2 for the dotted, 

V 

solid, and dashed curves respectively. zq=2 fm and £ A = 1.86. 
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11. Sea-quark and gluon distributions of MRS-1 (solid curves) and KMRS-BO (dotted 
curves) . 

12. Dependence on input distributions (see Fig. 11). Results by using the KMRS-BO 
inputs for (a) 12 C and (b) 40 Ca. Q 2 =0.8, 5, and 20 GeV 2 for the dotted, solid, and 
dashed curves respectively. zq=2 fm. £^=1.60 for C and 1.86 for Ca. <3 2 =0.8 GeV 2 . 
Compare them with the results in Figs. 7b and 8b. 

13. Dependence on the cutoff zq. Only Q 2 —5 GeV 2 curves are shown. zo=2.0 (3.0) fm 
for the solid (dashed) curves. C A = 1.86 and Q 2 =0.8 GeV 2 . 

14. Rescaling are used for all partons. Notations and parameters are same in Fig. 7b. 
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